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STEPS  INTO  COMPUTATIONAL  GEOMETRY: 
NOTEBOOK  II 


Abstract 

In  this  notebook  we  present  a collection  of  new  results  In 
computational  geometry,  which  all  concern  problems  of  planar  geometry. 

The  first  problem  Is  that  of  triangulating  a simple  n-vertex  polygon; 
we  show  that  this  can  be  done  In  time  O(nlogn),  by  first  decoiiq>oslng  In 
time  O(nlogn)  the  given  polygon  Into  a collection  of  special  polygons, 
called  monotone,  which  can  be  Individually  triangulated  In  time  propor- 
tional to  their  ntmibers  of  edges.  The  second  result  concerns  the  all- 
nearest  neighbor  problem  for  an  n-vertex  polygon:  a surprising  result 
Is  that,  also  no  method  faster  than  O(nlogn)  Is  known  for  constructing  the 
Voronol  diagram  of  a convex  polygon,  the  all-nearest-nelghbor  problem  can 
be  solved  In  time  0(n).  Finally,  we  show  the  feasibility  of  an  optimal 
real-time  algorithm  for  constructing  the  convex  hull  of  a set  of  n points 
In  the  plane.  This  algorithm  constructs  the  hull  by  successive  updates, 
using  total  O(nlogn)  time  - which  Is  optimal  - with  an  Interpolnt  delay 
O(logn),  which  gives  the  real-time  property. 


STEPS  INTO  COMPUTATIONAL  GEOMETRY 
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NOTEBOOK  II 

F.  P.  Preparata,  Editor 

Our  earlier  "Notebook",  by  the  title  of  "Steps  into  computational 
geometry",  was  issued  in  March  of  this  year  and  consisted  of  an  anthology 
of  selected  results  into  this  thriving  area  of  computational  complexity. 

After  a few  months,  we  have  assembled  a new  collection  of  results  to 
be  enclosed  in  this  report.  As  was  the  case  with  the  first  notebook,  there 
is  no  strong  unifying  scheme  for  the  results  to  be  presented,  with  the 
exception  that  they  all  concern  problems  of  planar  geometry.  One  of  the 
main  reasons  of  this  collection  is  simplicity  of  access  by  interested 
readers . 

The  first  problem  is  that  of  triangulating  a simple  n-vertex  polygon; 
we  show  that  this  can  be  done  in  time  0(nlogn},  by  first  decomposing  in 
time  O(nlogn)  the  given  polygon  into  a collection  of  special  polygons, 
called  monotone,  which  can  be  individually  triangulated  in  time  proportional 
to  their  numbers  of  edges.  The  second  result  concerns  the  all-nearest 
neighbor  problem  for  an  n-vertex  polygon:  a surprising  result  is  that,  also 
no  method  faster  than  O(nlogn)  is  known  for  constructing  the  Voronoi 
diagram  of  a convex  polygon,  the  all-nearest-neighbor  problem  can  be  solved 
in  time  0(n).  Finally,  we  show  the  feasibility  of  an  optimal  real-time 
algorithm  for  constructing  the  convex  hull  of  a set  of  n points  in  the 
plane.  This  algorithm  constructs  the  hull  by  successive  updates,  using 
total  O(nlogn)  time  - which  is  optimal  - with  an  interpoint  delay  O(logn), 
which  gives  the  real-time  property. 


TRIANGULATING  A SIMPLE  POLYGON 


Franco  P,  Prepara ta 

September  27,  1977 

The  efficient  algorithmic  construction  of  a trlangulatlon  of  a set  of 
points  In  the  plane  Is  an  Interesting  geometric  problem,  which  has  received 
considerable  Importance  from  recent  development  In  finite  element  methods  and 
Interpolation  techniques  11,2].  In  fact,  if  a function  f of  two  variables 
X and  y has  been  evaluated  at  a finite  set  S of  points  and  an  approximation 
Is  desired  at  a new  point.  It  Is  convenient  to  visualize  the  diagram  of  f 
as  a surface  consisting  of  triangular  plane  facets.  Therefore,  given  a 
trlangulatlon  of  S,  the  function  f can  be  evaluated  by  linear  Interpolation. 

The  problem  of  triangulating  a set  S of  n points  has  been  elegantly  solved 
by  Shamos  [3,41.  Indeed,  a trlangulatlon  Is  the  d\2al  graph  of  the  Voronol 
diagram  of  S [5],  which  is  a well-known  construct  for  solving  proximity 
problems  and  can  be  algorithmically  constructed  with  O(nlogn)  steps  on  a 

random-access  machine  with  real  number  arithmetic  [3,4]. 

2 

A more  difficult  problem  - for  which  only  0(n  ) brute  force  solutions  are 
known  so  far  - Is  the  trlangulatlon  of  a simple  polygon,  which  we  state  as 
follows:  "Given  an  n-vertex  simple  polygon  P,  subdivide  Its  Interior  Into 
triangles  whose  vertices  are  also  vertices  of  the  polygon". 

Notice  that  the  general  method  developed  for  a set  of  points  [3]  Is  not 
applicable  to  this  problem,  not  only  because  no  trlangulatlon  edges  may  exist 


This  work  was  supported  In  part  by  the  National  Science  Foundation  under 
Grant  MCS76-17321  and  in  part  by  the  Joint  Services  Electronics  Program 
under  Contract  DAAB-07-72-C-0259. 
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la  the  exterior  of  P,  but  also  because  the  edges  of  F must  belong  to 
the  trlangulatlon.  In  this  note  we  show  that  the  problem  can  also  be 
solved  with  at  most  O(nlogn)  steps. 

The  solution  will  be  obtained  In  two  stages,  which  are  kept  separate 
for  expository  reasons,  although  they  are  algorlthmi  ^lly  comblnable.  lu 
the  first  stage  we  decompose  In  time  at  most  O(nlogn)  the  given  simple 
polygon  P Into  a collection  of  simple  polygons  Pj^.P^, . . . , P^,  of  a special 
type  called  monotone.  Since,  as  we  shall  show,  a monotone  simple  polygon 
can  be  triangulated  In  time  proportional  to  the  number  of  its  vertices.  In 
the  second  step  we  triangulate  each  of  the  P^,  thereby  obtaining  the 
desired  trlangulatlon  of  P. 

For  convenience,  and  without  loss  of  generality,  we  choose  the  y-axls 
as  a preferred  direction.  We  say  that  a polygonal  line  or  chain  whose 
vertices  are  the  sequence  (Uj^,U2, . . . ,Up)  is  monotone  (with  respect  to  the 
y-axis)  if  y(Uj)  ^ y(u2)  ^ ...  y(Up).  A simple  polygon  P is  monotone  if  its 
boundary  consists  of  two  monotone  chains  with  common  extremes.  We  will  now 
prove  the  following  proposition. 

Proposition.  In  n-vertex  monotone  polygon  P can  be  triangulated  in  time 
0(a),  and  this  Is  optimal. 

Proof;  The  proof  Is  algorithmic.  The  algorithm  consists  of  at  most 
n-2-fv  steps,  where  v Is  the  number  of  nonconvex  (reflex)  vertices  of  the 
polygon.  Each  step  runs  In  time  bounded  by  a constant  and  Is  either  a 
trlangulatlon  step,  where  a trlangulatlon  edge  is  created,  or  a data 
manipulation  step.  For  easy  reference,  we  distinguish  In  P a left  and 
a right  chain,  with  obvious  meanings. 
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(fe  assume  Inductively  that  at  the  i-th  step  the  algorithm  examines  a 

structure  consisting  of  two  monotone  chains  (u«,u^ ,u_, . . , ,u  ) and  (u-jW) 

with  s s 1,  y(u  ) ^ y(w)  and.  If  s ^ 2,  (u  ..u.u^  ,)  > tt  (1)  for  j»l,...,s-l. 
® j“^x  3 3*1 

Without  loss  of  generality  we  shall  assume  that  ) is  on  the  left  of 

(Uq,w)  (see  figure  1).  The  data  structure  which  realizes  (Uq,...,u^)  Is  a 

stack,  stored  In  an  array  A,  with  two  pointers  V and  L,  such  that  L « V+1 

and  a[l]  >■  u , a[v]  » u , . Let  £ be  the  successor  of  u In  the  original  left 
s s-l  s 

chain  of  the  polygon  P.  We  distinguish  the  following  cases: 

(1)  y(f)  > y(w)  and  ^ (Xu  u .)  > tt  (figure  la).  This  Is  not  a 
trlangulatlon  step.  Vertex  i Is  added  to  the  left  chain,  l.e.,  V — V+1, 

L •-  L+1,  and  a[l]  ♦-  £ (In  other  words,  £ Is  "pushed"  Into  the  stack). 

(2)  y(X)  > y(w)  and  ^ (Xu  u .)  < tt  (figure  lb).  This  Is  a trlangula- 

tlon  step.  In  fact,  triangle  Xu  u , cannot  contain  any  vertex  of  the  polygon 

since  the  chains  are  monotone.  Thus,  the  trlangulatlon  edge  Xu  can  be 

S * X 

created.  At  this  point  we  check  ^(Xu  ^u  ,);  If  this  angle  is  > tt,  then 
we  update  the  stack  by  setting  aCl]  •"  X;  otherwise  we  set  V »-  V-1,  and 
create  the  trlangulatlon  edge  Xu  This  process  is  repeated  until  either 

<5  (XUj^u^  j^)  > tt  or  aCv]  =*  Uq  (the  stack  is  empty);  at  this  point  we  set 
L *-  V+l>  and  then  A[l]  •“  X. 

(3)  y(X)  ^ y(w)  (figure  Ic).  In  this  case  the  region  enclosed  by  the 
polygon  u.wu  u ^...u^  does  not  contain  any  vertex  of  the  polygon.  Moreover, 

by  the  hypothesis  that  <5  ^ ^ vertices 

u.,...,u  are  visible  from  w.  So  trlangulatlon  edges  wu, , for  j=l,...,s 

IS  j 

can  be  created. 


(a,b,c)"  denotes  the  counterclockwise  angle  formed  by  segments  ba 


and  be 
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Clearly,  q-3  trlangulatloa  edges  are  created  In  this  process,  each 
requiring  a fixed  amount  of  time;  case  (1)  can  occur  at  most  as  many  times 
as  there  are  reflex  verteces  in  P,  i.e.,  < n-3  times.  Thus  we  conclude 
that  P can  be  triangulated  in  time  0(n),  and  this  is  optimal. 


Notice  that  if  the  i-th  steps  examines  a two  chain  structure  as  described, 

(1)»  (2)»  (3)  above  so  will  the  (i+l)-st  step.  □ 

We  shall  now  show  that  any  simple  n-vertex  polygon  Q can  be  decomposed  in 
time  O(nlogn)  into  a collection  of  monotone  polygons,  by  means  of  a rather 
straightforward  modification  of  an  algorithm  by  Lee  and  Prepara ta  called 
"Regularization  of  a simple  polygon"  ([6],  p.603).  The  objective  of  the  latter 
procedure  was  the  addition  of  edges  to  a given  simple  polygon  to  ensure  that 
each  vertex  was  contained  in  some  monotone  chain  between  the  two  vertices  of 
largest  and  smallest  y-coordinate  (extreme  vertices).  This  was  done  as  follows. 
Assume,  for  simplicity  of  explanation,  that  no  two  vertices  in  the  vertex 

sequence  (Vq,Vj^ v^  j^)  of  the  polygon  Q have  identical  y-coordinates  (these 

ties,  howevei;  can  be  taken  care  of  very  easily).  Let  a cusp  be  a vertex  Vj 
for  which  either  y(Vj  j^)  < y(Vj)  and  y(Vj^j^)  < y(Vj)  or  y(Vj_j^)  > y(Vj)  and 
y(Vj^j^)  > y(Vj).  The  above  objective  is  met  if,  for  each  vertex  v of  Q - 
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except  the  two  extremes  vertices  - there  are  two  vertices  v'  and  v"  such 
that  y(v')  > y(v)  > y(  v") . This  condition  is  si ' Isf ied  by  all  vertices  except 
for  cusps:  thus  only  a cusp  - if  it  is  not  one  of  the  two  extreme  vertex  - 
is  to  be  connected  to  another  vertex  of  Q by  means  of  an  edge  not  crossing  any 
other  edge.  This  task  can  be  carried  out  by  the  above  mentioned  "regulariza- 
tion procedure"  in  time  O(logn)  per  vertex.  The  variant,  in  our 
present  problem  is  that  not  all  nonextreme  cusps  need  be  connected,  but 
only  those  whose  reflex  angle  is  inside  Q.  Thus  all  that  is  needed  is  a 
preliminary  scan  of  the  vertex  sequence  of  Q to  tag  only  the  cusp  which 
need  additional  connection;  next  we  can  apply  a regularization  procedure 
which  will  omit  connecting  untagged  vertices.  The  different  results  of  the 
original  [6]  and  of  the  modified  procedures  are  shown  in  Figures  2a  and  2b, 
respectively.  The  preliminary  scan  runs  in  time  0(n),  whereas  the  decompo- 
sition step  runs  in  time  O(nlogn).  The  latter  task  is  responsible  for  the 
order  of  complexity  of  the  entire  process.  In  figure  2(c)  we  show  the  end 
product  of  the  triangulation  process. 


Figure  2.  Simple  polygon  q:  (a)  regularized;  (b)  decomposed  into  monotone 
polygons;  (c)  triangulated. 
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THE  ALL  NEAREST-NEIGHBOR  PRDBIZM  FOR  CONVEX  POLYGONS* 

D.  T.  Lee 
August  1977 


Introduction 

The  problem  of  finding  the  nearest  neighbor  for  each  of  N arbitrary 
points  in  the  Euclidean  plane  has  been  shown  [l]  to  require  time  O(NlogN) , 
and  the  algorithms  achieving  the  lower  bound  have  also  been  given  in 
Cl,2j  . However,  the  lower  bound  does  not  apply  to  the  same  problem  when 
the  given  points,  rather  than  being  arbitrarily  placed,  form  a convex 
polygon.  In  this  paper,  we  shall  show  that  this  additional  information 
indeed  enables  us  to  obtain  a linear  time  algorithm,  whose  running  time  is 
obviously  optimal  within  a multiplicative  constant. 

Main  result 

Let  a convex  polygon  P be  denoted  by  a sequence  of  vertices  (PqjPj^ ^N-1^ 

in  which  p^p^^j, 0 * i < N,  is  an  edge.  Define  an  index  set  I={0, 1, . . . ,N-l} . 
Let  d(p^,p^),  i,j  € I,  denote  the  distance  between  p^  and  p^  and  D(P)  denote 
the  diameter  of  P,  i.e.,  D(P)  - ^ ^ ^ d(p^,Pj),  the  largest  distance  between 

This  work  was  supported  in  part  by  the  National  Science  Foundati on  under 
Grant  MCS-76-17321  and  in  part  by  the  Joint  Services  Electronics  Program 
under  Contract  DAAB-07-72-C-0259. 

^^Vll  indices  in  the  text  are  taken  modulo  N. 
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the  vertices  of  P.  The  nearest  neighbor  NN(p^)  of  is  p^  such  that 
d(Pj^>Pj)  “ g I-{i}*^^^i’^k^  * Consider  now  the  following  conditions; 

Condition  (i):  the  two  farthest  points  of  P are  the  extremes  of  an  edge, 

i.e.,  D(P)  = some  i. 

Condition  (ii):  all  vertices  of  P lie  inside  a circle  with  diameter  D(P) . 

A convex  polygon  P,  which  satisfies  both  (i)  and  (ii) , is  said  to  have  the 
semi-circle  property.  Figure  1 shows  a convex  polygon  having  the  semi-circle 
property. 


Figure  1.  A convex  polygon  with  the  semi -circle  property. 

Lemma  1 [3]  Given  a convex  polygon  P = (p^.p^^, . . . ,Pj^  , there  exists  a 

linear  time  algorithm  to  decompose  it  into  at  most  four  convex  polygons  which 
have  the  semi -circle  porperty. 

Proof:  First  of  all,  we  apply  the  linear  time  algorithm  [4]  to  find  the 
diameter.  Let  D(P)  = d(p^,p^).  The  chord  P^P^  will,  in  general,  divide  P 

into  two  convex  polygons  P^  - (Pu« Pu4-1’ * ' ' ’ V " 'Pn-I'^O ^u^ 

(Figure  2),  where  D(Pj^)  = D(P2)  = d(p^,p^).  Let  p^  € Pj^  be  the  vertex  with 
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largest  distance  from  the  chord  P^P^-  Let  be  defined  similarly. 

It  Is  obvious  that  p,  and  p can  be  found  In  0(N)  time.  p.  will  divide 

X n X 

Into  two  convex  polygons  • • -.P^)  ^12  Pf’PjJ+l' * ' ’'Pv^ 

Similarly,  p^  divides  P^  into  P^j  * (Py’^v+l’ * ‘ 

^22  * ^Pm’Piiri-l’"*’PN-l’P0’'**’Pu^*  claim  that  P^j^,  and  P^^ 

satisfy  the  semi-circle  property.  Without  loss  of  generality  we  shall  just 


fit 


Figure  2.  Decomposition  of  a convex  polygon  Into  four  convex  polygons 
satisfying  semi -circle  property. 


In  Figure  3,  since  p p is  the  longest  chord,  all  vertices  p ,P^,  ,..-,p 
must  lie  within  the  region  P^QP^  where  Q is  the  Intersection  of  the  two 
circles  with  radius  d(p^,p^)  and  centered  at  p^  and  p^  respectively.  By 
convexity  the  vertices  of  Pj^2  roust  lie  in  the  region  p^Q'p^,  and 


D(P--)  * d(p,,p  ).  The  line  A'p  perpendicular  to  p p Intersects  p.Q'  at  A'. 

IX  XV  V U V k 

The  region  p.Q'p  is  contained  in  the  right  triangle  p A'p.,  which  1 obviously 

At  yf  V X 

contained  In  the  semi -circle  with  diameter  p.p  and  center  A.  □ 

f V 


Figure  3.  Proof  that  P-,  has  semi -circle  property. 


Lemma  2 [3] ; Given  a convex  polygon  P = (p^.Pj^, . . . ,Pj^_j^)  with  the  semi-circle 
property,  for  any  vertex  p^^,  its  nearest  neighbor  Pj  is  adjacent  to  p^, 
i.e.,  either  j * i+1  or  j = i-1. 

Proof:  Without  loss  of  generality,  we  may  assume  that  D(P)  = • 

Suppose  for  some  p^^,  NN(pj^)  = where  k > i+1  (Figure  4)  . Consider  the 
triangle  Pj^Pj^^lPj^  • Since  ^Cp^^.P^)  ^^Pf^i+l^  assumption,  the  angle 

^ PiPi+lPk  angle  5 P^Pj^Pj^+l-  By  convexity,  p^  and  Pj^ 

must  lie  above  chords  PoP^+i  ^i+l^N-l  ‘ Th“®  ^ ^i^i+l^k 

must  be  greater  than  ^ PoPi+lPw-l  greater  than  rr/l  by  the  semi- 

circle property  of  the  given  polygon.  Ihat  is,  5 PiPfePi+l  ^ ^ ^i^i+l^k  ^ 
which  is  impossible.  Therefore  NN(Pj^)  must  be  adjacent  to  p^,  for 
all  i € I.  □ 
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% 

Figure  4.  Illustration  of  the  proof  of  lemma  2. 

Leiana  3.  Given  a convex  polygon  P • ^Po’^l ^N-1^  satisfying  con- 

dition (i)  above,  the  nearest  neighbor  NN(p^)  of  p^,  for  all  i € I,  can  be 
found  in  0(N)  time. 

Proof:  Suppose  D(P)  * dCp^,p„  .).  Let  p.  be  the  vertex  with  the  largest 
' " ~ U N-l  1 

distance  to  the  chord  By  lemma  1,  the  two  convex  polygons 

?!  = (Pq’Pi*  •••»  Pi^  ^2  “ ^^i’^i+l’  ^N-l^  satisfy  the  semi -circle 

property.  The  nearest  neighbor  of  each  vertex  in  P , (s  ■ 1,2)  can  be 

s 

fotind  separately  by  a simple  scan  through  the  vertices  of  P^  by  lemma  2. 

This  step  takes  0(N)  time.  We  must  still  check  whether  the  nearest  neighbor 
of  a vertex  in,  say  P^,  belongs  to  P^ , and  this  can  be  done  as  follows.  Let 
Pjj_j  be  the  origin,  and  the  chord  P^.j^Pq*  directed  from  p^_^  to  p^,  define 
the  positive  x-axis.  By  assumption,  p^  has  the  largest  y-coordinate . We 
project  all  the  vertices  on  the  vertical  line  p^q  through  p^  (Figure  5). 

The  projections  of  the  vertices  in  P^  and  P2,  denoted  by  * {j^CPq), 
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^(Pl),  •••.  ^(P^)}  and  /(P2)  “ {^(p^),  ^(P^+i).  •••.  are  ordered  as 

(■^(P^)y  ■*(Pj^_j)>  •••»  ^(PQ^y  (■*(Pj^)»  -^CPi+j) respectively, 

from  top  to  bottom.  Let 


Figure  5.  Projections  of  vertices  onto  pj^q  are  ordered  in  y-coordinate . 


6(Pj)  be  the  distance  from  Pj  to  its  nearest  neighbor  NN(Pj)  where  NN(Pj) 

and  p.  are  in  the  same  polygon  P_,  s = 1,2,  and  let  6-clrcle  (p.)  denote 
j a J 

the  circle  with  radius  6(p^)  and  centered  at  p^ . Consider  the  case  when 
Pj  € P^.  The  only  possible  nearest  neighbor  of  p^  is  among  those  Pj^  ^ ^2 
whose  projections  -^CPj^)  are  contained  in  the  S-circle  (p^)  (Figure  5). 

For  each  i(p^),  there  can  be  at  most  four  6-clrcles  that  pass  through  it 
[2].  Therefore,  in  order  to  find  NN(Pj)  for  all  Pj  ^ we  need  only 
examine  at  most  four  times  the  number  of  points  in  P2>  Specifically,  let 
the  vertices  of  P^  and  P2  be  kept  as  two  ordered  lists  (P^’P^.j*  •••*  Pq^ 

and  (Pj^»Pj^+l %-l^*  respectively  and  let  s^^  and  S2  be  pointers 

associated  with  these  two  lists,  respectively.  In  a first  scan  we  shall 
determine  the  nearest  neighbor  of  each  Pj€p^;  anologously,  in  a second 
scan  (not  to  be  described)  we  shall  process  P2*  In  each  step,  we  shall 


13 


examine  the  6-circle  (p„  ) for  p € p,  and  the  projection  i(p  ) for  p € P. . 

®1  ®1  ^ ®2  ®2  ^ 

We  start  by  setting  s,  ^ 1 and  s»  i+1.  If  the  6-circle  (p  ) does  not 

®1 

intersect  p.q  (Figure  5),  the  nearest  neighbor  of  p belongs  to  P.  and  was 
X S ^ X 

found  earlier.  So  we  update  s^  ^ s^  - 1.  Suppose,  for  some  p^  , the 

6-clrcle  (p  ) Intersects  p.q  at  u and  v where  u is  above  v (Figure  6). 

®1  ^ 

Then  we  examine  the  projection  X(p  ) to  see  if  It  belongs  to  the  segment 
— 2 

uv.  If  f (p  ) Is  above  u,  we  update  a.  ■«-  s.  + 1,  since  p can  not  be  the 
*2  ^ i S2 

nearest  neighbor  of  p . Let  f be  the  smallest  value  of  s„  so  that  /(p,) 

_ ^ ^ 
belongs  uv.  Since  p.  is  a possible  candidate,  we  check  if  d(p  ,p.)  is  less 

I S ^ X 

than  d(p  ,NN(p  )).  If  so,  we  update  NN(p  ) 

®1  ®1  ®1 
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8-  ♦ 8-  - 1.  Suppose  6-circle  (p  ,),  for  some  p , ^ P, , is  the  next  circle 

X X S ^ S ^ X 

that  intersects  p^^q  at  u'  and  v'  where  u'  is  above  v'  (Figure  6).  If  -6(p^) 

is  above  u' , we  advance  s^.  Otherwise,  we  scan  backward  and  if  -^Cp^)  is 

on  u'v*  we  do  the  checking  and  updating  as  before  until  we  reach  a vertex 

p whose  projection  i(p  ) is  above  u' ; then  we  scan  forward  from  f(p^)  on 

and  do  the  same  until  we  meet  a vertex  p^  whose  projection  ^^Cp^)  is  below 

v'.  At  this  point,  the  nearest  neighbor  of  p , is  determined,  we  update  s.. 

1 ^ 

Note  that  now  82  * r.  The  process  is  repeated  until  all  the  nearest 
neighbor  of  p^'s  in  are  determined.  As  pointed  out  above,  each  projec- 
tion can  only  be  examined  at  most  four  times,  so  the  time  required  for  the 
first  pass  is  0(n2)  where  n2  is  the  number  of  vertices  in  P2.  Similarly, 
the  time  for  the  second  pass  in  O(n^)  where  n^  is  the  number  of  vertices  in 
Pj^.  Therefore,  this  step  takes  0(nj+n2)  = 0(N)  time. 

Since  each  step  takes  0(N)  time,  the  proof  of  this  lemma  is 
completed.  □ 

Based  on  the  above  lemmas,  we  have  the  following  theorem. 

Theorem;  Given  a convex  polygon  P ■ (pQ,Pj,  ...»  Pjj.j)*  nearest 
neighbor  of  each  vertex  can  be  found  in  0(N)  time. 

Proof;  Let  D(P)  “ d(p^,p^).  The  chord  p^  p^  divides  the  polygon  P into 

two  polygons  P^  - (Pu»Pu+i*  Pv^  ^2  ’ ^^v’^v+l Pn-I’^O ^u^  • 

By  lemma  3,  the  nearest  neighbor  NN(p.)€P  of  p,  €P  , (s*l,2)  can  be 

J 8 j S 

fotmd  in  0(N)  time.  Now,  we  project  all  the  vertices  in  Pg,(s“l,2)  onto 
the  chord  P^P^'  Since  the  projections  of  the  vertices  in  P^,(8*l,2)  are 
ordered  respectively,  by  a technique  similar  to  that  descirlbed  in  lenna  3, 
we  can  find  for  each  vertex  p^  € P,  its  nearest  neighbor  in  0(N)  time. 


15 


Since  the  diameter  of  P can  be  found  in  0(N)  time,  the  total  running  time 
is  0(N).  □ 

Conclusion 

It  is  rather  interesting  that  the  nearest -neighbor  problem  for  a set 
of  N arbitrary  points  is  lower-bounded  by  0(N  log  N),  whereas  the  problem 
can  be  solved  in  linear  time  if  the  given  set  of  points  forms  a convex 
polygon.  In  [I],  the  nearest  neighbor  problem  was  solved  by  the  Voronoi 
diagram  technique.  The  construction  of  the  Voronoi  diagram  for  a set  of 
N points  has  also  been  shown  to  require  0(N  log  N)  time  [1].  But  whether 
the  construction  of  the  Voronoi  diagram  for  the  set  of  vertices  of  a 
convex  polygon  can  be  solved  in  less  than  0(N  log  N)  time  still  remains  an 
open  problem.  However,  we  know  at  least  that  the  nearest  neighbor  problem 
for  the  set  of  vertices  of  a convex  polygon  is  not  as  time-consuming  as 
the  presently  known  techniques  for  constructing  the  Voronoi  diagram  for  it. 
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AN  OPTIMAL  REAL-TIME  ALGORITHM  FOR  PLANAR  CONVEX  HULLS 
F.  P.  Preparata 


October  3,  1977 


1.  Introduction 

Algorithms  for  finding  the  convex  hull  of  a finite  set  of  n 
points  in  the  plane  have  been  developed  by  several  authors  in  recent 
years  [1,2, 3, 4].  ^fost  of  these  algorithms  are  also  optimal,  that  is,  as 
pointed  out  in  [5],  they  have  worst-case  running  time  0(n  log  n),  which  is 
also  the  best  achievable  performance. 

A common  feature  of  the  above  mentioned  algorithms  is  that  they 
are  all  off-line,  i.e.,  they  operate  on  the  data  collectively.  In  other 
words,  information  about  all  points  of  the  set  must  be  available  before 
any  of  those  algorithms  can  be  applied. 

Instead,  it  is  desirable  to  develop  an  algorithm  which  receives  one 
point  at  a time  and  updates  the  convex  hull  accordingly,  so  that,  after 

points  pj^,  p^ Pj^  have  been  received  their  convex  hull  is  available. 

Such  an  algorithm  is  appropriately  called  on-line.  A general  feature  of 
on-line  algorithms  is  that  no  bound  is  placed  on  the  update  time,  or, 
equivalently,  a new  item  (point)  is  input  on  request  as  soon  as  the  update 
relative  to  the  last  item  has  been  completed.  We  shall  refer  to  the  time 


This  work  was  supported  in  part  by  the  National  Science  Foundation  under 
Grant  MCS76-17321  and  in  part  by  the  Joint  Services  Electronics  Program 
under  Contract  DAAB-07-72-C-0259 . 
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Interval  between  two  consecutive  Inputs  as  the  interarrival  delay. 

Frequently,  known  on-line  algorithms  are  less  efficient  on  the 
entire  set  than  the  corresponding  off-line  algorithms  (some  price  must 
generally  be  paid  to  acquire  the  on-line  property).  For  the  planar  convex 
hull  problem,  however,  Shamos  has  designed  an  elegant  on-line  algorithm  [6], 
which  runs  In  time  O(nlogn),  thereby  matching  the  performance  of  off-line 
algorithms  for  the  same  problem. 

A more  demanding  case  of  on-line  applications  occurs  when  the  Inter- 
arrival delay  Is  outside  the  control  of  the  algorithm.  In  this  case  the 
update  must  be  completed  In  time  no  greater  than  the  minimum  Interarrlval 
delay.  Algorithms  for  such  applications  are  appropriately  called  in  real- 
time. Shamos  points  out  In  [6]  that,  since  any  convex  hull  algorithm  on 
n points  requires  H (nlogn)  operations,  any  real-time  algorithm  for  this 
problem  must  be  allowed  O(logn)  processing  time  between  successive  Inputs. 

Unfortunately,  the  algorithm  described  by  Shamos  exceeds  this 

2 

allowance,  since  its  Interarrlval  delay  can  be  0((logn)  ) . The  algorithm 
works  as  follows.  When  the  point  p^^  is  supplied,  assume  inductively,  that 
the  algorithm  has  available  the  convex  hull  of  the  set  of  points 

[p^,  ...,  a point  0 Internal  to  and  the  polar  angles  of  the 

vertices  of  " a convex  polygon  - about  0.  The  vertices  of  are 

arranged  In  a height-balanced  tree  (e.g.,  an  AVL  tree),  in  the  order 
of  their  polar  angles.  Thus  point  p^^  can  be  located  between  two  consecutive 
vertices  of  in  time  at  most  0(logi)  and  tested  for  Inclusion  in 

If  it  is  internal,  it  is  discarded;  otherwise,  two  vertices  I,  and  r of 
have  to  be  located  so  that  the  segments  p^jK  and  p^r  belong  to  lines  of 
support  of  y The  points  ^ and  r can  each  be  located  by  performing  a 
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standard  binary  search  on  the  vertex  cycle  of  other  hand,  each 

probe  of  this  search  is  Itself  a search  In  the  AVL  tree,  thereby  resulting 

2 2 

In  a worst-case  running  time  O((logi)  ) » 0((logn)  ) for  an  update. 

Ihe  Intuitive  reason  why  the  algorithm  sketched  above  falls  to 
achieve  an  O(logn)  update  time  is  that  a binary  search  Is  artificially 
forced  on  a search  tree,  rather  than  letting  the  latter  be  the  guide  of  the 
search  operation.  This  tiatural  observation  is  the  basis  of  the  following 
convex  hull  algorithm,  which  runs  in  time  O(nlogn),  and  is  therefore 
optimal,  and  has  update  time  0(logn},  and  is  therefore  in  real-time. 

2 . The  Real-time  Algorithm 

Let  P be  a polygon  in  the  plane  and  let  (Vq ^s-1^ 

counterclockwise  cycle  of  its  vertices  (Indices  are  modulo  s).  Ihe  vertices 
of  P will  be  stored  in  a data. structure  T(P)  which  is  a height  balanced  tree 
modified  in  a trivial  way.  Specifically  in  the  node  associated  with  vertex 
v^  we  also  store  a pointer  NEXT  [v^J , which  gives  the  address  of  the  node  of 
v^^j.  Let  also  min  T(P)  denote  the  first  member  of  the  vertex  cycle.  The 
convex-hull  algorithm  will  make  use  of  two  procedures:  TEST  and  RESTRUCTURE. 
Procedure  TEST  (P,m,p}  accepts  as  its  inputs  a point  p,  a convex  polygon  P, 
represented  by  the  tree  T(P),  and  the  minimum  element  m = min  T(P),  this 
algorithm  tests  whether  p is  internal  or  external  to  P,  and,  in  the  latter 
case,  it  determines  two  vertices  i and  r,  previously  defined  (see  Figure  1): 
notice  that  Jt  and  r are  named  so  that  ^ (rpj^)  < If  p is  internal,  the 

algorithm  terminates  without  altering  T(P);  otherwise  the  string  of  vertices 
comprised  between  I and  r is  deleted,  and  the  vertex  p is  then  inserted 
between  I and  r.  This  operation  is  performed  by  procedure  RESTRUCTURE  (P,p,i,r) 

(abc)  denotes  the  counterclockwise  angle  formed  by  segments  ba  and  be. 


Figure  1.  Defialeion  of  vertices  I and  r. 

Less  informally,  we  have  (A  is  the  empty  symbol): 

CONVEX-HULL  UPDATE 

Input;  ’^(^i-l^ » Pi 

Output;  T(Hj^) 

1.  begin  m *-  min 

2.  (Jt,r)  - TEST  (\_^,m.Pj^) 

3.  If  (jt,r)  ^ (A,A)  then  - RESTRUCTURE  j^.P^.i.r)  else 
end 

Obviously,  Step  1 runs  in  time  at  most  O(logl)  (search  in  height  balanced  tree 
with  at  most  i elements).  We  shall  now  show  that  both  TEST  and  RESTRUCTURE 
run  in  time  at  most  O(logl).  We  begin  by  considering  TEST.  Let  T « T(H^^j^), 
m ■ min(T),  and  M ■ ROOT(T).  Given  point  p^^  and  a vertex  v of  shall 

say  that  v is  convex  (with  respect  to  p^)  if  the  segment  p^v  intersects  the 
interior  of  otherwise,  if  the  two  vertices  adjacent  to  v lie  on  the 

same  side  of  the  line  containing  p^v,  v is  supporting;  in  the  re -mining  case, 

V is  reflex  (see  figure  2).  We  also  denote  as  a the  angle  (mp^M);  obviously 


Figure  2.  Illustration  of  convex,  supporting,  and  reflex  vertices. 


a la  classifiable  as  convex  rr)  or  reflex  (>  tt)  . Depending  upon  the 
classifications  of  m,  M,  and  at,  we  have  In  total  18  possible  cases.  These 
elementary  cases  can  be  conveniently  grouped  Into  eight  new  cases.  Illustrated 
In  figure  3,  each  requiring  distinct  algorithmic  actions.  Notatlonally,  let 
the  two  supporting  vertices  i and  r be  on  the  left  and  on  the  rlg^t  of  an 
observer  placed  in  and  facing  the  polygon  We  shall  make  use  of  two 

procedures,  called  LEFTSEARCH  and  R16HTSEARCH,  which  determine  I and  r, 
respectively.  The  arguments  of  these  procedures  are  binary  search  trees. 

We  also  let  L(M)  R(M)  denote  the  left  and  the  right  subtree  of  the  node  M, 
respectively.  The  analysis  of  the  eight  cases  and  their  corresponding  actions 
Is  straightforward.  It  Is  worth  pointing  out  that  when  p^  is  Internal  to 
cases  1 or  7 will  occur  repeatedly,  the  algorithm  will  examine  a nested  family 
of  subtrees, and  will  terminate  when  the  subtree  consists  of  only  one  leaf, 
i.e.,  when  m - M (see  Step  2 below).  We  can  now  explicitly  give  the  algorithm: 

TEST  (H,m,p) 

Input:  H,  a polygon,  represented  as  a modified  AVL  tree  T(H} 
p a point,  m the  minimum  element  In  T(H). 

Output : either  a pair  (i,r)  of  Integers  or  (A, A) 

1.  begin  M - ROOT  (T(H)),  T *-  T(H) 

2.  ^ m ■ M then  v *-  t •-  fi  (Cocmiant;  p^  Is  Internal) 

3.  else  begin  ^ a ^ n then 

4.  ^ m Is  convex  then 

5.  ^ M is  convex  then  T - R(M),  m - NEXT(M),  u - 0 (case  1) 

6.  else  Tj  - T-R(M) , T2  •-  R(M) , u •-  1 (case  2) 

7.  else  ^ M is  reflex  then  T *-  L(M) , u •-  0 (case  3) 

8.  else  Tj^  •“  T - L(M) , T^  *"  I*(M) , u *-  I (case  4) 

9.  else  If  m Is  reflex  then 

^ M Is  reflex  then  T •-  R(M),m  •-  NEXT(M)  u •-  0 (case  5) 
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11. 

12. 

13. 

14. 

15. 


16. 


and 

return  (i,r) 


else  Tj^  *-  R(M),T2  - T - R(M) , u *-  1 (case  6) 
else  ^ M is  convex  then  T - L(M) , u •-  0 (case  7) 
else  - L(M),  T2  - T - L(M) , u 1 (case  8) 
If  u » 0 then  (f,r)  *-  TEST(T,in,Pj^) 
else  I - LEFTSEARCH(Tj^),  r - RIGHTSEARCH(T2) 


end 


L(M) 


R(M) 


Figure  3.  The  eight  possible  cases  handled  by  algorithm  TEST. 


We  shall  now  describe  the  procedure  LEFTSEARCH  used  by  the  algorithm 
TEST  (RIGHTSEARCH  is  analogous  with  obvious  modifications) . 

Procedure  LEFTSEARCH 


Input;  a tree  T,  describing  a sequence  of  vertices 
Output ; a vertex  I 

1.  begin  c *-  ROOT(T) 

2.  ^ pc  is  supporting  then  I — c 

3.  else  begin  ^ c is  reflex  then  T L(c)  else  T R(c) 

4.  I - LEFTSEARCH(T) 
end 

5.  return  I 


-and- 


! 
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It  is  obvioxis  that  LEFTSEARCH  involves  tracing  a path  of  the  tree  T,  j 

spending  a bounded  time  at  each  node.  Since  T is  a balanced  tree  with  at 

most  (i-1)  nodes,  the  mnning  time  is  O(logi).  Referring  now  to  the  time 
performance  of  TEST,  we  notice  that  the  bulk  of  the  work  is  done  either  in 
Step  14  (recursive  call  of  TEST)  or  in  Step  15  (calls  of  LEFTSEARCH  AND 
RIGHTSEARCH) , whereas  the  decisions  leading  to  either  of  these  steps  (Steps 
3-13)  take  time  bounded  by 'a  constant.  Typically,  algorithm  TEST  could  be  viewed  as 
tracing  a path  from  the  root  to  some  node  c of  T(H^  > recursively  calling 
itself.  If  p^  is  internal  to  then  c is  a leaf  of  T(H^  j^)  ; otherwise, 

starting  at  node  c,  two  paths  of  T(H^_j^)  are  traced  by  LEFTSEARCH  and 
RIGHTSEARCH,  respectively,  until  X and  r are  found.  Since  the  amount  of  work 
expended  at  each  node  is  bounded  by  a constant,  TEST  runs  in  time  O(logi). 

Finally,  we  consider  the  procedure  RESTRUCTURE,  which  is  invoked  only 
when  p^  is  external  to  H^  Let  n^^  ^ be  the  number  of  vertices  ^ . | 

As  mentioned  earlier,  the  vertices  comprised  between  X and  r must  be  deleted 
and  p^  must  be  inserted.  With  regard  to  the  deletion,  slightly  different  actions  | 
will  be  taken  depending  upon  whether  X precedes  r in  T(H^  ^)  or  not.  In  the 
first  case  we  have  to  split  twice  and  splice  once  AVL  trees  with  at  most  I 

i-1  elements;  in  the  second  case,  only  two  splittings  occur.  But  split  and  | 

splice  of  AVL  trees  are  standard  operations,  known  as  Crane's  algorithms  ^ 

([7],  p.465),  which  can  be  performed  in  time  O(logi)  and  will  not  be  further  j 

discussed.  Similarly, insertion  of  p^  can  be  done  in  time  O(logi),  whereas  the 
update  of  the  function  NEXT  only  Involves  two  pointers,  associated  with  X ^ 

and  p^  respectively. 

Therefore,  we  conclude  that  CONVEX -HULL  UPDATE  can  be  executed  in  time 
O(logi)  after  i points  have  been  processed  and  can  be  used  for  an  optimal 
real-time  algorithm. 
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